Environmental  Laboratory  erdc/el tr-16-ii 


US  Army  Corps 
of  Engineers® 

Engineer  Research  and 
Development  Center 


INNOVATIVE  SOLUTIONS 
for  a  safer,  better  world 


Ecosystem  Management  and  Restoration  Research  Program 

Testing  and  Validation  Studies  of  the  NSMII- 
Benthic  Sediment  Diagenesis  Module 


Zhonglong  Zhang  and  Billy  E.  Johnson 


July  2016 


J3 

C 

QJ 

E 


Algae ;  POC, 

DO  demand 

PON f,  POP ;  Si 

CH4 

h2s  r 

Deposition 

i 

A 

V 

NH4 


N03 


DIP 


DSi 


Transfer 


CH4 


H2S(d) 


NH4  (d) 


N03 


POCqi,2 


POC 


POCG2,2 


Djf  fusion 

Diagenesis 


cm 


H2S(d) 


pocg1 


pong1i2 


PON 


PONG2/2 


Diagenesis 


NH4  (d) 


PONG3,2 


N03 


P04  (d) 


T 


Si(d) 


POPG1,2 


POP  I POP G2,2  _ 


Diagenesis 

w 

▼ 

P04  (d) 

W 

POPG3/2 


Partitioning 


H2S(p) 


Particle 
mixing 
Partitioning 


article  t 
nixing  j 


OqO 
OpO 
CH4,  H2S,  yQO 

NH3,  CO 2  ° 

OqO 

o  o 

ttt 

OqO 
O  O 

o°R°o 

o°o 

o§°8 
— §30- 

o2°o 


NH4  (p) 


P04  (p) 


Si  (p) 


h2s(p) 


Partitioning 


NH4  (p) 


T 


Partitioning 


P04  (p) 


T 


BSi 


Dissolution 


Burial 


Si(d) 


Partitioning 


Si  (p) 


~T~ 

Burial 


Approved  for  public  release;  distribution  is  unlimited. 


The  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC)  solves 
the  nation’s  toughest  engineering  and  environmental  challenges.  ERDC  develops 
innovative  solutions  in  civil  and  military  engineering,  geospatial  sciences,  water 
resources,  and  environmental  sciences  for  the  Army,  the  Department  of  Defense, 
civilian  agencies,  and  our  nation’s  public  good.  Find  out  more  at 
www.erdc.usace.army.mil. 


To  search  for  other  technical  reports  published  by  ERDC,  visit  the  ERDC  online  library 
at  http://acwc.sdp.sirsi.net/client/default. 


Ecosystem  Management  and 
Restoration  Research  Program 


ERDC/EL  TR- 16-11 
July  2016 


Testing  and  Validation  Studies  of  the  NSMII 
Benthic  Sediment  Diagenesis  Module 


Zhonglong  Zhang 

LimnoTech 
501  Davis  Drive 
Ann  Arbor,  Ml  48108 

Billy  E.  Johnson 
Environmental  Laboratory 

U.S.  Army  Engineer  Research  and  Development  Center 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 


Final  report 

Approved  for  public  release;  distribution  is  unlimited. 


Prepared  for  U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 

under  Project  45745,  “Nutrient  Simulation  Module  Development  and  Integration" 


ERDC/EL  TR-16-11 


ii 


Abstract 

Nutrient  simulation  modules  (NSMs)  have  been  developed  as  “plug  in” 
water  quality  modules  for  HEC-RAS  (Hydrologic  Engineering  Center- 
River  Analysis  System)  and  other  hydrologic  and  hydraulic  models.  The 
NSMs  model  multiple  water  quality  constituents  and  biogeochemical  pro¬ 
cesses  in  the  water  column.  The  NSMs  consist  of  two  modules:  NSMI  and 
NSMII,  differing  mainly  in  the  number  of  water  quality  constituents  (state 
variables)  and  kinetic  processes  (transformations)  included.  The  benthic 
sediment  diagenesis  module  was  developed  as  an  option  module  for  cou¬ 
pling  with  the  NSMI  and  NSMII  water  column  kinetics.  Sediment-water 
fluxes  of  dissolved  oxygen  and  nutrients  are  computed  internally  rather 
than  prescribed  if  the  benthic  sediment  diagenesis  module  is  activated. 

The  ability  of  the  NSMII-benthic  sediment  diagenesis  module  to  correctly 
predict  sediment-water  nutrient  fluxes  and  sediment  oxygen  demand  was 
validated  and  evaluated  against  observed  data  and  existing  benthic  sedi¬ 
ment  diagenesis  models  through  a  series  of  case  studies.  These  case  stud¬ 
ies  were  chosen  for  representing  various  sedimentation  and  environmental 
conditions  in  aquatic  systems.  The  simulated  sediment-water  nutrient 
fluxes  and  sediment  oxygen  demand  over  time  were  generally  in  good 
agreement  with  observed  data  and  two  model  results  for  all  test  cases.  This 
report  documents  testing  and  evaluation  of  the  NSMII-benthic  sediment 
diagenesis  module.  The  report  describes  the  testing  conditions,  model  in¬ 
puts,  and  model  results. 
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1  Introduction 

1.1  Background 

A  set  of  nutrient  simulation  modules  (NSMs)  have  been  developed  as  “plug 
in”  water  quality  modules  for  HEC-RAS  (Hydrologic  Engineering  Center- 
River  Analysis  System)  (HEC  2010)  and  other  hydrologic  and  hydraulic 
models  (Zhang  and  Johnson  2016).  The  NSMs  are  designed  to  be  inde¬ 
pendent  of  the  dimensionality  of  the  spatial  domain;  they  model  multiple 
water  quality  constituents  and  biogeochemical  processes  in  aquatic  sys¬ 
tems.  The  NSMs  consist  of  two  kinetics  modules:  NSMI  and  NSMII,  differ¬ 
ing  mainly  in  the  number  of  water  quality  constituents  (state  variables) 
and  kinetic  processes  (transformations)  included.  They  are  compiled  as 
NSMI.dll  and  NSMII.dll,  respectively.  The  NSMII  was  designed  for  an  ad¬ 
vanced  water  quality  simulation.  The  NSMII  models  multiple  algal  groups, 
benthic  algae,  nitrogen,  phosphorus,  and  carbon  cycles,  carbonaceous  bio¬ 
chemical  oxygen  demand,  dissolved  oxygen,  and  pathogen.  The  constitu¬ 
ents  simulated  in  NSMII  are  listed  in  Table  1.  Chla,  D,  C,  N,  P  and  O2 
under  the  units  refer  to  chlorophyll-a,  dry  weight,  carbon,  nitrogen,  phos¬ 
phorus,  and  oxygen,  respectively.  The  colony  forming  unit  (CFU)  is  a 
measure  of  viable  bacterial  numbers. 


Table  1.  Water  quality  state  variables  modeled  in  NSMII. 


Variable 

Definition 

Units 

Option 

A  pi 

Algae  (Phytoplankton) 

pg-Chla  L1 

1-3 

At, 

Benthic  algae 

g-D  nr2 

On/Off 

N03 

Nitrate 

mg-N  L1 

On 

NH4 

Ammonium 

mg-N  L1 

On 

DON 

Dissolved  organic  nitrogen 

mg-N  L1 

On 

LPON 

Labile  particulate  organic  nitrogen 

mg-N  L1 

On 

RPON 

Refractory  particulate  organic 
nitrogen 

mg-N  L1 

On 

TIP 

Total  inorganic  phosphorus 

mg-P  L 1 

On 

DOP 

Dissolved  organic  phosphorus 

mg-P  L 1 

On 

LPOP 

Labile  particulate  organic 
phosphorus 

mg-P  L 1 

On 

RPOP 

Refractory  particulate  organic 
phosphorus 

mg-P  L 1 

On 

DIC 

Dissolved  inorganic  carbon 

mol  L1 

On 
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Variable 

Definition 

Units 

Option 

LDOC 

Labile  dissolved  organic  carbon 

mg-C  L1 

On 

RDOC 

Refractory  dissolved  organic  carbon 

mg-C  L 1 

On 

LPOC 

Labile  particulate  organic  carbon 

mg-C  L1 

On 

RPOC 

Refractory  particulate  organic 
carbon 

mg-C  L1 

On 

CBOD, 

Carbonaceous  biochemical  oxygen 
demand 

mg-02  L1 

0-10 

CH4 

Methane 

mg-02  L1 

On/Off 

HxS 

Total  dissolved  sulfides 

mg-02  L1 

On/Off 

DO 

Dissolved  oxygen 

mg-02  L1 

On 

BSi 

Particulate  biogenic  silica 

mg-Si  L 1 

On/Off 

DSi 

Dissolved  silica 

mg-Si  L1 

On/Off 

PX 

Pathogen 

cfu  (100  mL)1 

On/Off 

Aik 

Alkalinity 

mg-CaC03  L1 

On/Off 

Figure  1  provides  an  overview  of  the  water  quality  state  variables  and  mass 
transport  processes  modeled  in  the  NSMII.  Algae  may  be  either  floating 
(phytoplankton)  or  attached  to  the  bed,  in  a  manner  similar  to  periphyton. 
Floating  algae  are  subject  to  sinking  while  periphyton  are  subject  to  sub¬ 
strate  limitations.  Diatoms  are  distinguished  from  other  algae  in  that  they 
need  silicate  to  grow.  Organic  matter  has  a  fundamental  role  in  water 
quality  processes,  which  cannot  be  identified  by  single  lumped  organic  ni¬ 
trogen,  phosphorus,  and  carbon  state  variables  alone  (Connolly  and  Coffin 
1995,  Shanahan  et  al.  1998,  Chapra  1999).  Compared  to  NSMI,  NSMII  in¬ 
cludes  dissolved,  refractory,  and  labile  particulate  organic  species.  Labile 
and  refractory  organic  matter  distinctions  are  based  upon  relative  decay 
rates  of  the  organics.  A  labile  fraction  describes  organic  materials  that  de¬ 
cay  on  a  time  scale  of  days  to  weeks,  while  a  refractory  fraction  accounts 
for  decay  processes  lasting  months  to  a  year.  Note  that  the  NSMII  splits 
dissolved  organic  carbon,  particulate  organic  carbon,  nitrogen,  and  phos¬ 
phorus  into  labile  and  refractory  fractions  but  represents  dissolved  organic 
nitrogen  and  phosphorus  as  homogenous  components.  Water  quality  ki¬ 
netics  for  CBOD,  pathogen,  and  alkalinity  are  modeled  in  NSMII  as  for 
NSMI.  Detailed  discussion  of  equations  and  mechanisms  in  NSMI  and 
NSMII  can  be  found  from  Zhang  and  Johnson  (2016). 
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Figure  1.  Diagram  of  water  quality  state  variables 
and  mass  transport  processes  in  NSMII. 


Sediment  oxygen  demand  (SOD),  due  to  the  diagenesis  of  organic  materi¬ 
als  in  bottom  sediments,  can  contribute  to  oxygen  declines  in  water  bod¬ 
ies,  such  as  seasonal  declines  in  hypolimnetic  oxygen  observed  in  most 
reservoirs  (Martin  et  al.  2007).  Over  lengthy  time  scales  (e.g.,  years  to  dec¬ 
ades),  the  benthic  sediments  are  an  ultimate  sink  of  nutrients  and  other 
substances  discharged  to  the  water  column.  Over  lesser  time  scales  (e.g., 
seasons  to  years),  however,  sediment  release  of  previously  deposited  nu¬ 
trients  can  be  a  net  source  to  the  water  column  (Di  Toro  et  al.  1990,  Chap- 
ra  et  al.  2015).  Modeling  SOD  and  nutrient  releases  at  the  sediment-water 
interface  is  fundamental  to  understanding  and  mitigating  the  eutrophica¬ 
tion  and  hypoxia  problems  in  water  bodies.  Therefore,  an  optional  benthic 
sediment  diagenesis  module  was  included  for  coupling  with  the  NSMII’s 
water  column  kinetics.  Sediment-water  fluxes  of  dissolved  oxygen  and  nu¬ 
trients  are  computed  internally  rather  than  prescribed  if  the  benthic  sedi¬ 
ment  diagenesis  module  is  activated. 

1.2  Objective 

A  thorough  validation  and  evaluation  of  the  newly  developed  sediment  di¬ 
agenesis  module  is  essential.  The  authors  were  not  aware  of  measurements 
that  allow  a  direct  validation  of  the  results  of  the  benthic  sediment  diagen- 
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esis  module.  Therefore,  the  sediment  diagenesis  module  developed  for 
NSMI  and  NSMII  has  been  verified  and  evaluated  through  a  series  of  case 
studies.  In  these  test  cases,  sediment  diagenesis  simulation  computed  re¬ 
sults  were  compared  to  similar  outputs  generated  from  the  sediment  dia¬ 
genesis  model  including  Microsoft  Excel  SedFlux  (Chapra  et  al.  2008),  and 
CE-QUAL-ICM  (Cerco  and  Cole  1993,  Di  Toro  and  Fitzpatrick  1993,  Cerco 
et  al.  2004).  This  report  describes  the  results  of  the  model  validation  and 
evaluation  studies.  The  report  also  provides  information  about  the  sedi¬ 
ment  diagenesis  input  parameters  and  model  capabilities. 

1.3  Approach 

In  this  study,  intercomparison  was  carried  out  between  three  different  wa¬ 
ter  quality  models,  each  with  a  sediment  diagenesis  component,  including 
Microsoft  Excel  SedFlux  (Chapra  et  al.  2008),  and  CE-QUAL-ICM  (Cerco 
and  Cole  1993,  Di  Toro  and  Fitzpatrick  1993,  Cerco  et  al.  2004).  The  inter¬ 
comparison  included  a  comparison  between  model  results  and  algorithms. 
A  series  of  case  studies  were  designed  for  representing  various  sedimenta¬ 
tion  and  environmental  conditions.  Note  that  validation  in  this  context  re¬ 
fers  to  validation  of  the  model  framework,  rather  than  validation  of 
simulation  results.  Model  validation  tests  were  designed  to  analyze  benthic 
sediment  diagenesis  processes  and  estimate  sediment-water  exchanges  of 
nutrients  and  oxygen  demand.  All  test  case  studies  were  used  to  verify  that 
the  algorithms  and  formulations  are  implemented  correctly  in  NSMII. 
Model  testing  and  validation  studies  were  also  used  to  provide  further  in¬ 
sight  into  the  sediment-water  interactions  between  the  processes  regulat¬ 
ing  the  SOD  and  nutrient  release.  This  report  documents  testing  and 
evaluation  of  the  NSMII-benthic  sediment  diagenesis  module. 
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2  Benthic  Sediment  Diagenesis  Module 

The  NSMs  couple  the  water  column  kinetics  with  a  two-layer  benthic  sed¬ 
iment  diagenesis  module  and  internally  computes  sediment-water  column 
fluxes  of  nutrients  and  SOD  based  on  deposition  of  organic  particles  and 
other  factors.  The  basic  framework  of  the  sediment  diagenesis  module 
consists  of  two  well-mixed  sediment  layers,  a  thin  upper  layer  (layer  l)  and 
a  thicker  active  layer  on  the  order  of  10  cm  (layer  2).  The  upper  layer  (layer 

1) ,  in  contact  with  the  water  column,  maybe  oxic  or  anoxic,  depending  on 
DO  concentration  in  the  water.  The  lower  layer  (layer  2)  is  always  anoxic. 
The  upper  layer  depth  is  at  its  maximum  only  a  small  fraction  of  the  active 
sediment  depth  (~  0.1  cm).  The  sediment  diagenesis  module  simulates 
four  basic  processes:  1)  the  deposition  of  algae  and  particulate  organic 
matter  (POC,  PON,  POP)  as  depositional  fluxes  from  the  water  column  di¬ 
rectly  to  the  second  layer  due  to  the  negligible  thickness  of  the  upper  layer, 

2)  the  diagenesis  of  the  sediment  particulate  organic  matter  in  the  second 
layer  produces  soluble  substances  that  are  quantified  as  diagenesis  fluxes, 

3)  the  solutes  react  and  transfer  between  two  layers  or  are  released  as  gas¬ 
es,  and  4)  the  solutes  are  returned  to  the  overlaying  water.  Figure  2  pro¬ 
vides  an  overview  of  the  state  variables  and  processes  simulated  in  the 
sediment  diagenesis  module.  Diagenesis  (mineralization)  reactions  occur 
in  the  second  layer.  Organic  matter  initially  mineralizes  rapidly  in  the  ben¬ 
thic  sediments,  but  then  slows  down.  In  order  to  capture  this  process,  sed¬ 
iment  POC,  PON,  and  POP  are  fractionated  to  three  “G  classes.”  The  basis 
of  the  multi-G  assumption  came  from  laboratory  experiments  that  showed 
organic  matter  decay  could  be  approximated  as  a  function  of  different 
pools  (Berner  1980).  The  Gi,  labile,  fraction  has  a  half-life  of  20  days.  The 
G2,  refractory,  fraction  has  a  half-life  of  one  year.  The  G3,  inert,  fraction 
undergoes  no  significant  decay  before  burial  into  deep,  inactive  sediments 
(Di  Toro  2001). 
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Figure  2.  Schematic  overview  of  the  state  variables  and  processes  modeled  in  the 
benthic  sediment  diagenesis  module. 


Additionally,  27  sediment  diagenesis  state  variables  are  included  in  each 
NSM  module  if  the  sediment  diagenesis  module  is  activated.  Table  2  lists 
the  sediment  diagenesis  state  variables  and  the  symbols.  Three  state  vari¬ 
ables  are  utilized  to  represent  three  reactivity  classes  (Gi  through  G3).  The 
units  of  these  variables  are  grams  per  sediment  layer  volume  basis,  or  mg 
Lr1.  Methane,  sulfate,  and  sulfide  are  tracked  in  units  of  oxygen  equivalents 
to  easily  balance  the  model’s  computations. 


Table  2.  Water  quality  state  variables  modeled  in  the  NSMIl-benthic  sediment 

diagenesis  module. 


Symbol 

Variable  Name 

Bed  layer 

Units 

NH4i 

Sediment  ammonium 

1 

mg-N  L1 

N03i 

Sediment  nitrate-nitrite 

1 

mg-N  L1 

CH4i 

Sediment  methane 

1 

mg-02  L1 

S04i 

Sediment  sulfate 

1 

mg-C>2  L 1 

TH2S1 

Sediment  total  sulfide 

1 

mg-02  L1 

TIPi 

Sediment  total  inorganic  phosphorous 

1 

mg-P  L1 

DSii 

Sediment  dissolved  silica 

1 

mg-Si  L1 
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Symbol 

Variable  Name 

Bed  layer 

Units 

NH42 

Sediment  ammonium 

2 

mg-N  L1 

N032 

Sediment  nitrate-nitrite 

2 

mg-N  L1 

CH42 

Sediment  methane 

2 

mg-02  L 1 

SO  4  2 

Sediment  sulfate 

2 

mg-02  L1 

TH2S2 

Sediment  total  sulfide 

2 

mg-02  L1 

tip2 

Sediment  total  inorganic  phosphorous 

2 

mg-P  L1 

DS/2 

Sediment  dissolved  silica 

2 

mg-Si  L1 

P0Car2 

Sediment  particulate  organic  carbon  (G1 
-  G3) 

2 

mg-C  LA 

PONgi.2 

Sediment  particulate  organic  nitrogen 
(G1  -  G3) 

2 

mg-N  L1 

POPg/,2 

Sediment  particulate  organic 
phosphorous  (G1  -  G3) 

2 

mg-P  L1 

BSh 

Sediment  particulate  biogenic  silica 

2 

mg-Si  L1 

ST 

Sediment  benthic  stress 

2 

d 

The  sediment  diagenesis  module  employed  here  includes  two  options  for 
its  numerical  solutions:  l)  steady  state  and  2)  unsteady  state.  Initial  condi¬ 
tions  are  required  for  the  unsteady-state  solution  for  each  of  the  sediment 
cells  below  the  water  column  cells.  Initial  conditions  are  specified  sepa¬ 
rately  for  the  sediment  cells,  and  differ  from  those  for  the  water  column. 
Initial  conditions  are  required  only  for  state  variables  at  sediment  layer  2 
listed  above.  A  detailed  discussion  of  the  theory  and  formulations  in  the 
benthic  sediment  diagenesis  module  can  be  found  in  Zhang  and  Johnson 
(2016). 
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3  Testing  and  Validation  Studies 

A  comprehensive  validation  of  the  benthic  sediment  diagenesis  module 
would  require  an  extensive  and  contemporary  field  measurement  of  all  the 
variables  considered.  There  would  be  many  difficulties  involved  with  col¬ 
lecting  the  necessary  information.  However,  some  independent  checks  on 
the  various  coefficients  and  formulations  are  desirable.  Thus,  a  set  of  mod¬ 
el  testing  cases  were  designed  to  investigate,  primarily,  the  model  perfor¬ 
mance.  The  sediment  diagenesis  module  was  calibrated  and  validated 
against  experimental  data  (Martin  2002)  and  the  Chesapeake  Bay  dataset 
from  1991  to  2000  through  a  series  of  case  studies.  Since  the  benthic  sed¬ 
iment  diagenesis  module  is  based  upon  the  same  framework  as  Microsoft 
Excel  SedFlux  (called  SedFlux  henceforth)  and  CE-QUAL-ICM  (called  ICM 
henceforth),  model  results  from  SedFlux  and  ICM  were  used  in  the  follow¬ 
ing  model  validation  studies.  The  SedFlux  is  a  standalone  model  of  simu¬ 
lating  SOD  and  sediment-water  fluxes  of  nutrients.  The  subroutine  used  in 
SedFlux  was  adapted  from  original  work  for  the  QUAL2K  model  (Chapra 
et  al.  2008).  SedFlux  simulates  the  production  and  oxidation  of  sulfide  in¬ 
stead  of  methane  during  sediment  diagenesis  in  brackish,  estuarine,  or 
marine  environments  where  sufficient  sulfate  is  present.  SedFlux  includes 
a  steady-state  mode  and  a  time  variable  (unsteady)  mode.  The  ICM’s  sed¬ 
iment  diagenesis  model  was  originally  developed  as  a  standalone  code  and 
then  incorporated  into  the  ICM  model  (Di  Toro  and  Fitzpatrick  1993)  for 
the  Chesapeake  Bay  application.  This  sediment  diagenesis  module  was  fur¬ 
ther  modified  to  update  methanogenesis-sulfate  reduction  reactions  (Cer- 
co  et  al.  2004).  There  are  major  differences  in  the  model  codes  relating  to 
the  computation  of  SOD  and  methane  for  the  three  models.  The  ICM  mod¬ 
el  includes  methane  and  sulfides  in  the  computation  of  SOD,  but  differs  in 
the  implementation  of  the  sulfide/methane  reactions.  In  SedFlux,  sulfide 
is  computed  for  the  saltwater  case,  while  methane  is  computed  for  the 
freshwater  case.  Sulfides,  sulfate,  and  methane  are  computed  in  the  ICM 
and  NSMII  for  both  conditions,  with  a  regression  equation  used  to  esti¬ 
mate  sulfates  based  upon  salinity.  In  SedFlux,  sediment  methane  concen¬ 
tration  is  only  computed  using  a  steady-state  solution  for  the  aerobic  layer. 
Silica  is  modeled  in  the  sediment  diagenesis  module,  but  is  not  included  in 
these  testing  and  validation  studies. 
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A  testing  and  validation  study  was  conducted  using  the  Chesapeake  Bay 
water  quality  datasets  through  comparing  simulation  results  generated  by 
the  NSMII,  with  model  results  generated  with  SedFlux  and  ICM.  The  three 
models  were  run  with,  to  the  extent  possible,  a  consistent  set  of  kinetic  co¬ 
efficients  and  parameters,  initial  conditions,  deposition  fluxes  from  the 
water  column  and  water  column  concentrations  in  all  test  cases.  In  partic¬ 
ular,  NSMII  and  SedFlux  sediment  diagenesis  parameters  and  coefficients 
for  the  model  comparison  were  adopted  from  a  calibrated  ICM  model  ap¬ 
plication  to  the  Chesapeake  Bay.  The  NSMII-benthic  sediment  diagenesis 
module  was  further  tested  and  validated  through  several  special  cases  and 
a  sensitivity  analysis. 

3.1  Model  development  and  result  comparisons 

In  this  test  case,  an  equivalent  water  column  and  benthic  sediment  seg¬ 
ment  was  set  up  for  SedFlux,  ICM  and  NSMII.  The  sediment  diagenesis 
computational  time  step  was  set  equal  to  the  water  quality  model  time 
step.  Model  predictions  by  the  NSMII  were  compared  against  model  re¬ 
sults  generated  from  SedFlux  and  ICM  using  the  1991-2000  Chesapeake 
Bay  dataset.  The  Chesapeake  Bay  dataset  was  used  to  test  and  validate  the 
performances  of  the  NSMII  predictions  with  respect  to  various  aspects 
such  as  sediment  layer  organic  and  inorganic  concentrations,  SOD,  and 
sediment-water  transfer  fluxes  of  inorganic  nutrients. 

3.1.1  Input  data  and  model  development 

A  computational  water  quality  cell  with  the  same  force  conditions  was 
used  to  develop  three  models.  Forcing  conditions  included  the  sediment 
deposition  rate,  a  determinant  of  sediment  organic  matter  supply  and  the 
sediment  concentration,  and  thus  the  driver  for  diagenesis.  Solids  concen¬ 
trations  in  sediment  layers  1  and  2  were  set  to  equal  0.5  kg  Lr1.  Thickness 
of  active  sediments  was  equal  to  10  cm.  Specified  initial  concentrations  of 
sediment  diagenesis  state  variables  are  given  in  Table  3. 

Table  3.  Initial  conditions  of  benthic  sediment  diagenesis  state  variables. 


Symbol 

Bed  layer 

Units 

Values 

NH42 

2 

mg-N  L1 

7.76 

N032 

2 

mg-N  L1 

0.07 

CH42 

2 

mg-02  L1 

0.02 

S042 

2 

mg-C>2  L1 

1031.67 
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Symbol 

Bed  layer 

Units 

Values 

TH2S2 

2 

mg-02  L1 

3590.40 

TIP2 

2 

mg-P  L1 

112.47 

P0Cgi,2 

2 

mg-C  L 1 

56.23 

P0Cg2,2 

2 

mg-C  L1 

1072.33 

P0Cg3,2 

2 

mg-C  L 1 

24966.55 

P0Ngi,2 

2 

mg-N  L1 

10.87 

P0Ng2,2 

2 

mg-N  L1 

208.34 

P0Ng3,2 

2 

mg-N  L1 

3114.92 

P0Pgi,2 

2 

mg-P  L1 

0.83 

POP  G2,2 

2 

mg-P  L 1 

9.25 

P0PG3,2 

2 

mg-P  L 1 

98.30 

ST 

2 

d 

22.55 

The  NSMII  computes  time-varying  concentrations  in  a  water  column  over- 
lying  the  sediment  bed.  In  order  to  identify  differences  in  the  time-variable 
nature  of  model  predictions  for  the  three  models,  The  Chesapeake  Bay  da¬ 
taset  used  in  the  ICM  model  inputs  were  read  into  both  the  NSMII  and 
SedFlux.  This  test  case  ensured  that  all  three  models  used  the  same  forcing 
functions  for  deriving  the  bed  sediment  diagenesis  processes.  Depositional 
rates  of  water  column  organic  matter  to  the  benthic  sediments  are  given  as 
model  inputs  and  shown  in  Figure  3. 

Figure  3.  Time  series  plots  of  water  column  depositional  rates: 

(a)  POC,  (b)  PON,  and  (c)  POP. 
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Water  column  concentrations  of  ammonia  (NH4),  nitrate  (NO3),  dissolved 
inorganic  phosphorous  (DIP),  and  dissolved  oxygen  (DO)  are  shown  in 
Figure  4.  DO  time  series  exhibit  summer  peaks  at  this  location  where  over- 
lying-water  DO  levels  remain  above  10  mg/L  for  the  majority  of  the  year. 


Figure  4.  Time  series  plots  of  the  water  column  concentrations: 
(a)  NH4,  (b)  N03,  (c)  DIP,  and  (d)  DO. 
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Unsteady-state  SedFlux,  NSMII  and  ICM  models  were  run  with  the  condi¬ 
tions  specified  above  and  results  compared.  Each  model  was  run  for  a  pe¬ 
riod  of  looo  days.  Table  4  presents  the  final  values  of  the  sediment 
diagenesis  process  parameters  and  coefficients  used  in  the  three  models. 
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Table  4.  Key  sediment  diagenesis  module  input  parameters  and  coefficients. 


Symbol 

Parameter 

Bed 

layer 

Unit 

Value 

Temperature 

dependent 

o 

o 

Q 

Layer  1  critical  O2  for  incremental  sorption 

1 

mg-02  L1 

2.0 

A/CP041 

Layer  1  increamental  inorganic  P  partition 
coefficient 

1 

unitless 

300 

Vch4,l(T) 

Layer  1  CH4  oxidation  transfer  velocity 

1 

m  d1 

0.2 

Yes 

1.079 

Vh2s,d(T) 

Layer  1  dissolved  H2S  oxidation  transfer 
velocity 

1 

m  day1 

0.2 

Yes 

1.079 

Vh2s,p(T) 

Layer  1  particulate  H2S  oxidation  transfer 
velocity 

1 

m  day1 

0.4 

Yes 

1.079 

ksH2S 

Layer  1  H2S  oxidation  normalization 
constant 

1 

mg-02  L 1 

4.0 

KsOxch 

Layer  1  half-saturation  oxygen  constant 
for  CH4  oxidation 

1 

mg-02  L1 

0.2 

Vnh4,l(T) 

Layer  1  nitrification  transfer  velocity 

1 

m  d1 

0.14 

Yes 

1.08 

KsOxnal 

Layer  1  half-saturation  oxygen  constant 
for  nitrification 

1 

mg-02  L 1 

0.37 

KsNh4 

Layer  1  half-saturation  NH4  constant  for 
nitrification 

1 

mg-N  L 1 

0.728 

Vno3,l(T) 

Layer  1  denitrification  transfer  velocity 

1 

m  d1 

0.125 

Yes 

1.08 

Css 

Layer  1  solid  concentration 

1 

kg  L1 

0.5 

Fapi 

Fraction  of  algae  settling  into  G1  class 

2 

unitless 

0.6 

FaP2 

Fraction  of  algae  settling  into  G2  class 

2 

unitless 

0.2 

Frpoci 

Fraction  of  RPOC  settling  into  sediment 
P0CG1 

2 

unitless 

0.15 

F RP0C2 

Fraction  of  RPOC  settling  into  sediment 
P0CG2 

2 

unitless 

0.35 

Frponi 

Fraction  of  RPON  settling  into  sediment 
PONG1 

2 

unitless 

0.15 

F RP0N2 

Fraction  of  RPON  settling  into  sediment 
PON  G2 

2 

unitless 

0.35 

Frpopi 

Fraction  of  RPOP  settling  into  sediment 
POP  G1 

2 

unitless 

0.15 

F RP0P2 

Fraction  of  RPOP  settling  into  sediment 
POP  G2 

2 

unitless 

0.35 

Kpocgi(T) 

Diagenesis  rate  of  sediment  POC  G1 

2 

d1 

0.035 

Yes 

1.1 

Kp0CG2(T) 

Diagnesis  rate  of  sediment  POC  G2 

2 

d1 

0.0018 

Yes 

1.15 

Kpongi(T) 

Diagenesis  rate  of  sediment  PON  G1 

2 

d1 

0.035 

Yes 

1.1 

KpONG2(T) 

Diagenesis  rate  of  sediment  PON  G2 

2 

d1 

0.0018 

Yes 

1.15 
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Symbol 

Parameter 

Bed 

layer 

Unit 

Value 

Temperature 

dependent 

Kpopgi(T) 

Diagenesis  rate  of  sediment  POP  G1 

2 

d1 

0.035 

Yes 

1.1 

Kp0PG2(T) 

Diagenesis  rate  of  sediment  POP  G2 

2 

d1 

0.0018 

Yes 

1.15 

kdnh42 

Layer  2  NH4  partition  coefficient 

2 

L  kg1 

1.0 

kdh2s2 

Layer  2  H2S  partition  coefficient 

2 

L  kg1 

100 

Vno3,2(T) 

Layer  2  denitrification  transfer  velocity 

2 

m  d1 

0.25 

Yes 

1.08 

kdpo42 

Layer  2  inorganic  P  partition  coefficient 

2 

L  kg1 

100 

kst 

Decay  rate  of  sediment  benthic  stress 

2 

d1 

0.03 

KsS04 

Layer  2  S04  half-saturation  constant  for 
reduction 

2 

mg-02  L1 

1.08 

Css 

Layer  2  solid  concentration 

2 

kg  L1 

0.5 

W2 

Sedimentation  burial  velocity 

2 

m  d1 

4.93E- 

05 

Dd(T) 

Sediment  pore-water  diffusion  coefficient 

m2  d1 

0.0005 

Yes 

1.08 

Dp(T) 

Sediment  particle  mixing  diffusion 
coefficient 

m2  d1 

0.00003 

Yes 

1.117 

KsDp 

Half-saturation  oxygen  constant  for 
sediment  particle  phase  mixing 

mg-02  L1 

4.0 

POCr 

Reference  sediment  POC  for  bioturbation 

mg-C  g1 

0.1 

3.1.2  Predictions  of  sediment  organic  matter  diagenesis  fluxes 

Sediment  organic  matter  diagenesis  processes  are  modeled  by  partitioning 
the  settling  organic  matter  into  three  reactivity  classes  (3G).  Each  class 
represents  a  fixed  portion  of  the  organic  material  that  reacts  at  a  specific 
rate.  Diagenesis  of  particulate  organic  carbon  (POC),  particulate  organic 
nitrogen  (PON),  and  phosphorous  (P)  is  computed  from  a  mass  balance 
equation.  Model  predictions  for  sediment  diagenesis  fluxes,  e.g.,  carbon 
(Jcc),  nitrogen  (Jn),  and  phosphorous  (Jp),  are  shown  in  Figure  5.  The 
three  models  (SedFlux,  ICM,  and  NSMII)  produced  identical  results  as  ex¬ 
pected  because  the  same  water  quality  deposition  inputs  and  diagenesis 
rates  were  applied  in  the  three  models. 
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Figure  5.  Comparisons  of  model  predictions  of  the  sediment  total  diagenesis  fluxes: 
(a)  carbon,  (b)  nitrogen,  (c)  phosphorous. 


(a) 


(b) 


(c) 


3.1.3  Predictions  of  sediment  reaction  and  transfer  coefficients 

Dissolved  (Kli2)  and  particle  (W12)  mixing  between  layers  1  and  2  are  modeled  as 
a  function  of  passive  transport  and  proxies  for  the  activities  of  benthic  organisms. 
The  rate  of  mixing  of  sediment  particles  (W12)  is  quantified  by  estimating  the  ap¬ 
parent  particle  diffusion  coefficient.  Dissolved  species  are  transported  through  the 
two  sediment  layers  by  diffusion  (Km).  The  sediment-water  mass  transfer  coeffi¬ 
cient  (Kloi)  controls  solute  exchange  between  the  aerobic  layer  and  the  overlying 
water  column.  All  three  parameters  depend  on  the  thickness  of  the  active  sedi¬ 
ment  layer  and  are  internally  calculated  in  the  NSMII. 
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Model  predictions  for  the  particle  mixing  coefficient  (W12),  diffusion  coefficient 
(KL12),  sediment-water  transfer  coefficient  (Kloi),  and  benthic  stress  (ST)  are 
shown  in  Figure  5.  The  three  models  (SedFlux,  ICM,  and  NSMII)  produced  the 
identical  Kli2  as  expected  because  the  same  formulation  is  implemented.  Values 
of  W12,  Kloi,  and  ST  computed  from  the  ICM  and  NSMII  are  approximately  the 
same.  The  modeled  values  varied  slightly  between  the  NSMII  and  SedFlux.  This 
is  because  the  sediment-water  transfer  and  benthic  stress  are  computed  by  the 
SedFlux  differently  from  the  other  two  models.  In  the  SedFlux,  either  methane  or 
sulfide  is  used  to  compute  carbonaceous  oxygen  demand  (CSOD),  which  is  less 
than  the  CSOD  computed  from  the  ICM  and  NSMII.  The  smaller  SOD  causes  a 
smaller  sediment-water  transfer  coefficient  shown  in  Figure  6(c).  The  temperature 
correction  on  the  benthic  stress  is  not  included  in  the  SedFlux. 


Figure  6.  Comparisons  of  model  predictions  of  (a)  particle  mixing  coefficient  (W12),  (b) 
diffusion  coefficient  (Kli2),  (c)  sediment-water  transfer  coefficient  (Kloi),  and  (d) 

sediment  stress  (ST). 


(a) 


(b) 
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3.1.4  Predictions  of  sediment  concentrations  of  inorganic  constituents 

Sediment  concentrations  of  ammonia  (NH4),  nitrate  (NO3),  total  inorgan¬ 
ic  phosphorous  (TIP),  sulfide  (TH2S),  and  methane  (CH4)  are  computed 
for  both  aerobic  and  anaerobic  layers  via  mass  balances  of  biogeochemical 
and  physical  processes.  Model  predictions  for  sediment  NH4,  NO3,  TIP, 
TH2S,  and  CH4  in  layers  2  and  1  are  shown  in  Figure  7  and  Figure  8,  re¬ 
spectively.  Because  the  same  sediment  diagenesis  fluxes  were  produced  as 
shown  in  Figure  5,  the  three  models  (SedFlux,  ICM,  and  NSMII)  produced 
approximately  the  same  sediment  concentrations  except  sulfide  and  me¬ 
thane.  In  the  SedFlux,  sulfide  is  computed  for  the  saltwater  case,  while 
methane  is  computed  for  the  freshwater  case.  In  the  ICM  and  NSMII,  sul¬ 
fate,  sulfide,  and  methane  are  computed  for  both  conditions.  Sediment 
sulfate  is  included  as  a  state  variable  in  the  ICM  and  NSMII.  Sulfate  is 
produced  in  the  aerobic  layer  as  a  result  of  particulate  and  dissolved  sul¬ 
fide  oxidation  and  removed  in  the  anaerobic  layer  via  sulfate  reduction. 
Sulfate  limitation  on  sulfide  production  (via  sulfate  reduction)  is  modeled. 
Methanogenesis  is  computed  as  a  fraction  of  total  carbon  diagenesis, 
which  is  also  controlled  by  anaerobic  layer  sulfate.  Methanogenesis  is 
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equal  to  total  carbon  diagenesis  when  anaerobic  layer  sulfate  is  extremely 
low  relative  to  the  half-saturation  constant.  Therefore,  sediment  sulfide 
and  methane  predicted  by  the  SedFlux  are  different  from  the  results  com¬ 
puted  by  the  other  two  models. 


Figure  7.  Comparisons  of  model  predictions  of  the  sediment  layer  2’s  concentrations: 
(a)  NH4,  (b)  N03,  (c)  TIP,  (d)  H2S,  and  (e)  CH4. 
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Sediment  concentrations  for  ammonia  (NH4),  nitrate  (NO3),  total  inor¬ 
ganic  phosphorous  (TIP),  sulfide  (TH2S),  and  methane  (CH4)  in  layer  1 
follow  the  same  pattern  as  layer  2.  The  three  models  (SedFlux,  ICM,  and 
NSMII)  produced  approximately  the  same  sediment  concentrations  except 
sulfide  and  methane.  These  differences  are  attributed  to  different  algo¬ 
rithms  implemented  by  each  model. 


Figure  8.  Comparisons  of  model  predictions  of  the  sediment  layer  l’s  concentrations: 
(a)  NH4,  (b)  N03,  (c)  TIP,  (d)  H2S,  and  (e)  CH4. 
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3.1.5  Predictions  of  sediment  oxygen  demand  and  sediment-water 
fluxes 

SOD  is  computed  by  summing  all  model  processes  that  consume  oxygen, 
including  the  oxidation  of  NH4,  H2S,  and  CH4.  The  CSOD  is  determined 
in  the  ICM  and  NSMII  from  both  methane  and  sulfide  oxidation.  In  the 
SedFlux,  either  methane  or  sulfide  is  used  to  compute  the  CSOD.  Sedi¬ 
ment-water  fluxes  of  dissolved  variables  (NH4,  NO3,  DIP,  H2S,  and  CH4) 
are  determined  by  their  concentration  differences  between  the  surface  sed¬ 
iment  layer  and  the  overlying  water.  The  flux  may  be  in  either  direction 
across  the  sediment-water  interface,  depending  on  the  concentration  gra¬ 
dient.  Positive  fluxes  are  from  benthic  sediment  to  water  column.  Negative 
fluxes  are  from  water  column  to  sediments.  Model  predictions  of  SOD  and 
sediment-water  fluxes  for  NH4,  NO3,  DIP,  CH4,  and  H2S  are  shown  in 
Figure  9.  Seasonal  cycles  of  SOD  at  this  location  where  the  overlying  water 
column  may  be  oxic  and  hypoxic  are  out  of  phase  with  one  another.  Strong 
oxygen  limitation  of  SOD  occurs,  which  is  illustrated  by  occasional,  short¬ 
lived  SOD  peaks  during  summer  when  DO  concentrations  are  slightly  ele¬ 
vated.  SOD  computed  from  the  SedFlux  varied  slightly  between  the  Sed¬ 
Flux  and  ICM  and  NSMII  because  only  CSOD  derived  by  sediment  sulfide 
is  computed  in  the  SedFlux.  For  the  same  reason,  sediment-water  fluxes  of 
methane  and  sulfide  computed  from  the  SedFlux  are  different  from  the 
other  two  models.  Otherwise,  the  three  models  (SedFlux,  ICM,  and 
NSMII)  produced  approximately  the  same  results  for  sediment-water  flux¬ 
es  of  ammonia,  and  dissolved  inorganic  phosphorus. 


Figure  9.  Comparisons  of  SedFlux  and  NSMII  unsteady-state  model  predictions  of  the 
sediment-water  fluxes:  (a)  SOD,  (b)  NH4,  (c)  N03,  (d)  DIP,  (e)  CH4  and  (f)  H2S. 
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In  summary,  the  NSMII-benthic  sediment  diagenesis  module  was  able  to 
predict  POC,  PON,  and  POP  deposition  fluxes  into  each  G  class  in  the  sed¬ 
iment  layer.  Inorganic  nutrient  concentrations  in  sediment  layers  l  and  2 
derived  by  the  sediment  diagenesis  processes  were  consistently  simulated 
against  model  results  produced  by  the  SedFlux  and  ICM.  The  NSMII  pre¬ 
dictions  of  SOD  and  sediment-water  fluxes  of  nutrients  matched  reasona¬ 
bly  well  with  the  model  results  obtained  from  the  SedFlux  and  ICM.  Some 
differences  among  the  three  models  as  seen  are  due  to  different  algorithms 
implemented  in  each  model.  The  model  validation  was  successful. 

3.2  Special  case  testing  and  sensitivity  analysis  of  NSMII 

This  section  describes  several  special  testing  and  validation  case  studies  of 
the  NSMII’s  sediment  diagenesis  module. 

3.2.1  NSMII  vs.  SedFlux  steady-state  solutions 

Both  SedFlux  and  NSMII  include  two  options  for  the  numerical  solutions 
of  mass  balance  equations  of  sediment  diagenesis  state  variables:  l) 
steady-state  and  2)  unsteady-state.  The  first  option  does  not  require  initial 
conditions  for  state  variables.  Initial  conditions  are  required  for  the  second 
option.  This  test  case  provides  comparisons  of  model  predictions  given  by 
these  two  models  under  steady-state  conditions.  Concentrations  of  all  sed¬ 
iment  diagenesis  state  variables  predicted  by  these  two  models  are  identi¬ 
cal  except  methane  and  sulfide  concentrations.  The  predicted  methane 
and  sulfide  concentrations  varied  slightly  between  the  two  models.  In  the 
SedFlux,  sulfide  is  computed  only  for  the  saltwater  case,  while  methane  is 
computed  for  the  freshwater  case.  Sediment  methane  is  assumed  in  the 
SedFlux  to  be  completely  mixed  for  both  layers.  The  methane  concentra¬ 
tion  is  then  only  computed  under  steady-state  solutions.  Steady-state 
model  predictions  for  SOD  and  sediment-water  fluxes  of  ammonia,  ni- 
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trate,  phosphate,  and  methane  are  shown  in  Figure  10.  The  two  models 
produced  comparable  predictions  for  SOD  and  nutrient  release  rates. 
Temporal  variations  of  SOD  between  the  SedFlux  and  NSMII  are  slightly 
different  due  to  sediment  methane  concentrations. 


Figure  10.  Comparisons  of  SedFlux  vs.  NSMII  steady-state  solutions:  (a)  SOD  and 
sediment-water  fluxes:  (b)  NH4,  (c)  N03,  (d)  CH4,  and  (e)  DIP. 
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3.2.2  NSMII  steady-state  vs.  unsteady-state  solutions 

In  this  validation  case,  two  NSMII  model  simulations  were  conducted  for 
the  steady-state  and  unsteady-state  conditions.  Initial  sediment  concen¬ 
trations  of  all  state  variables  were  set  to  zero  for  one  simulation.  Model 
predictions  for  sediment  layer  2’s  ammonia  (NH4),  nitrate  (NO3),  dis¬ 
solved  inorganic  phosphorous  (DIP),  sulfide  (H2S),  and  methane  (CH4) 
are  shown  in  Figure  11.  The  results  produced  from  both  steady  and  un- 
steady-state  solutions  are  different.  Initial  sediment  concentrations  for 
NH4,  NO3,  DIP,  CH4,  and  H2S  under  unsteady-state  solutions  have  effects 
on  the  sediment-water  fluxes  at  the  beginning  of  the  simulation  period; 
however,  their  concentrations  reach  the  same  after  that.  Therefore,  initial 
sediment  concentrations  of  each  state  variable  from  unsteady-state  solu¬ 
tions  have  minimum  impact  on  the  long-term  concentration  of  each  state 
variable  predicted  by  the  NSMII. 

The  steady-state  mode  of  the  sediment  diagenesis  module  is  typically  used 
to  estimate  initial  conditions.  A  steady-state  computation  allows  the  model 
to  compute  the  initial  conditions  based  upon  the  initial  water  column  con¬ 
centrations  and  organic  matter  settling  fluxes,  assuming  steady-state  con- 
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ditions.  As  advised  by  Di  Toro  (2001),  a  common  approach  to  refine  the 
initial  conditions  is  through  the  use  of  a  quasi-dynamic  approach.  In  this 
approach,  the  model  is  run  for  a  period  of  a  year  or  more,  with  reasonably 
repeating  water  column  conditions,  and  the  predicted  concentrations  in 
diagenesis  segments  at  the  end  of  that  period  are  used  to  refine  and  re¬ 
place  the  specified  initial  conditions.  This  process  is  repeated  until  the  re¬ 
sulting  predictions  approach  a  quasi-steady  state. 


Figure  11.  Comparisons  of  NSMII  steady-  vs.  unsteady-state  solutions  of  sediment 
layer  2’s  concentrations:  (a)  NH4,  (b)  N03,  (3)  DIP,  (4)  H2S,  and  (c)  CH4. 
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3.2.3  NSMII  analytical  vs.  numerical  solutions  of  sediment  methane 

The  dissolved  methane  produced  by  carbon  diagenesis  can  be  oxidized  in 
the  aerobic  layer,  released  as  a  flux  to  the  overlying  water  column,  or  be¬ 
come  gaseous  when  its  concentration  exceeds  saturation.  The  methane 
saturation  is  used  to  establish  an  upper  limit  for  the  sediment  methane 
concentration  in  layer  2  and  corresponds  to  the  CSOD  that  results  from 
the  methane  oxidation.  If  the  methane  resulting  from  the  sediment  carbon 
diagenesis  is  greater  than  the  methane  produced  at  saturation,  the  CSOD 
is  computed  from  the  saturation  value  and  the  remaining  methane  may  be 
stored  and  removed  as  gas  bubbles  (Rudd  and  Taylor  1980).  The  loss  of 
methane  as  bubbles  is  an  important  sink.  For  this  validation  test,  the  POC 
depositional  flux  from  the  water  column  was  increased  by  a  factor  of  100 
over  that  used  in  the  above  test  case.  This  test  was  included  to  allow  suffi¬ 
cient  POC  flux  to  the  sediment  bed  so  that  the  resulting  methane  produc¬ 
tion  would  exceed  the  saturated  concentration.  In  this  test  case,  modeled 
sediment  methane  concentration  in  layer  2  (CH42)  under  the  numerical 
solution  is  set  to  saturation  (CH4s)  as  shown  in  Figure  12(a).  In  the  NSMII, 
it  is  assumed  that  super-saturation  does  not  occur  and  that  all  methane 
produced  in  excess  of  the  saturation  concentration  is  immediately  trans- 
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ferred  to  gas  bubbles.  This  means  that  all  methane  produced  after  the  es¬ 
tablishment  of  the  saturation  concentration  is  lost  to  the  atmosphere. 


An  analytical  solution  is  used  to  determine  the  steady-state  flux  of  dis¬ 
solved  methane  corrected  for  gas  loss  delivered  to  the  aerobic  sediments 
and  the  oxygen  consumed  by  oxidation  of  methane.  In  this  approach,  sed¬ 
iment  methane  is  only  computed  in  layer  l.  Comparisons  of  the  gas  losses 
of  sediment  methane  as  bubbles,  CSOD,  and  sediment-water  flux  of  me¬ 
thane  under  analytical  and  numerical  solutions  are  shown  in  Figure  12(b), 
12(c),  and  12(d),  respectively.  CSOD  computed  from  the  numerical  solu¬ 
tion  is  lower  than  the  analytical  solution,  having  a  2  to  3  g-02/m3  differ¬ 
ence  at  the  CSOD  peak.  The  sediment-water  flux  of  methane  is  also  lower 
than  the  analytical  solution.  The  gas  losses  of  sediment  methane  computed 
from  the  numerical  solution  is  higher  than  the  analytical  solution.  If  the 
methane  saturation  is  exceeded,  the  numerical  solution  is  not  appropriate 
since  the  CSOD  becomes  constant  and  hits  a  ceiling  (Chapra  1997). 


Figure  12.  Comparisons  of  NSMII  analytical  vs.  numerical  solutions  of  sediment 
methane:  (a)  layer  2’s  CH4,  (b)  CSOD,  (c)  CH4  loss,  and  (d)  sediment-water  flux  of 
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3.2.4  Comparisons  of  the  diagenesis  rates  of  three  sediment  organic 
matter  classes  (Gl,  G2,  G3) 

Sediment  diagenesis  (mineralization)  reactions  occur  in  the  second  layer. 
The  heterogeneous  nature  of  sediment  organic  matter  is  recognized  and 
simulated  by  adopting  a  multi-G  approach  described  by  Di  Toro  (2001). 
Three  “G  classes”  (labile,  refractory,  inert)  are  modeled  in  the  NSMII.  La¬ 
bile,  refractory,  and  inert  distinctions  are  based  upon  the  time-scales  of 
oxidation  or  decomposition.  The  Gl  fraction  has  a  half-life  of  20  days.  The 
G2  fraction  has  a  half-life  of  one  year.  The  G3  fraction  undergoes  no  signif¬ 
icant  decay  before  burial  into  deep,  inactive  sediments.  In  this  validation 
case,  digenesis  process  of  the  three  G  classes  of  sediment  POC,  PON  and 
POP  is  tested  and  compared  with  a  steady-state  solution.  The  steady-  and 
unsteady-state  model  predictions  for  sediment  POC,  PON  and  POP  Gl  and 
G2  classes  are  shown  in  Figure  13.  The  POC,  PON  and  POP  G  class  1  ap¬ 
proached  the  steady-state  concentrations  relatively  quickly,  while  the  pre¬ 
dictions  for  G  class  2  did  not  reach  a  steady  state  by  the  end  of  the  1000- 
day  simulation.  This  test  case  indicates  that  it  is  necessary  to  divide  the 
sediment  organic  matter  into  multiple  groups. 
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Figure  13.  Comparisons  of  NSMII  predictions  of  the  sediment  concentrations  under 
the  steady-  and  unsteady-state  solutions:  (a)  POC  Gl,  (b)  POC  G2,  (c)  PON  Gl,  (d) 
PON  G2,  (e)  POP  Gl,  and  (f)  POP  G2. 
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3.2.5  Sensitivities  of  the  water  column  depositional  rates 

Sensitivity  analysis  is  an  important  part  of  the  model  testing  and  valida¬ 
tion.  In  this  validation  case,  sensitivity  analysis  of  depositional  rates  of  or¬ 
ganic  matter  from  the  water  column  was  used  for  determining  the  most 
important  forcing  functions  as  well  as  the  potential  influence  of  SOD  and 
release  of  phosphate,  ammonia,  and  nitrogen  from  the  sediment.  The  in¬ 
put  of  organic  matter  to  the  sediment-water  interface  is  one  of  the  most 
important  chemical  inputs  driving  the  sediment  diagenesis.  In  this  sensi¬ 
tivity  analysis,  water  column  concentrations  of  POC,  PON,  and  POP  speci- 


ERDC/EL  TR-16-11 


31 


fied  in  the  above  tests  were  reduced  by  20%,  50%,  and  90%,  respectively. 
Model  predictions  for  SOD  and  sediment  release  of  NH4,  NO3,  H2S,  CH4, 
and  DIP  are  shown  in  Figure  14.  The  results  show  that  reduced  deposi- 
tional  loads  produced  lower  SOD  and  sediment-water  fluxes  of  NH4,  NO3, 
H2S,  CH4,  and  DIP.  The  sediment  concentrations  of  both  organic  and  in¬ 
organic  substances  were  reduced  proportionally  to  depositional  fluxes 
from  the  water  column.  Therefore,  depositional  fluxes  of  organic  matter 
have  important  influence  on  the  benthic  sediment  concentrations  and  nu¬ 
trient  releases.  The  SOD  and  sediment  release  of  nutrients  are  most  sensi¬ 
tive  to  different  depositional  rates  without  considering  dynamic 
sedimentation  processes. 


Figure  14.  Sensitivities  of  water  column  organic  matter  deposition  rates  on  sediment- 
water  fluxes:  (a)  SOD,  (b)  NH4,  (c)  N03,  (d)  H2S,  and  (e)  DIP. 
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3.2.6  Sensitivities  of  the  deep  burial  velocities 

In  the  sediment  diagenesis  module,  the  active  sediment  layer  does  not 
move  vertically,  its  position  relative  to  the  sediment-water  interface 
changes  as  additional  material  is  deposited  on  the  bottom.  Thus,  deposi¬ 
tion  of  particulate  matter  at  the  sediment  water  interface  results  in  older 
sediments  moving  further  from  that  boundary.  The  resuspension  process 
results  in  older  sediments  moving  from  the  sediment-water  interface  and 
is  not  included  in  the  model.  Only  burial  mechanism  establishes  the  verti- 
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cal  position  of  constituents  in  the  sediment  profile.  In  applying  the  NSMII, 
the  burial  velocities  specified  in  terms  of  input  parameters  are  not  associ¬ 
ated  with  the  flux  of  total  suspended  solids  (inorganic  and  organic)  from 
the  overlying  water.  In  this  sensitivity  analysis,  the  burial  velocity  was  in¬ 
creased  by  a  factor  of  10  and  too  than  that  used  in  the  above  test  case  (w2 
=  4.928-10-5  m/d).  No  burial  was  also  included  to  test  the  sensitivity  of 
burial  velocities  on  the  SOD  and  sediment  release  of  nutrients.  Model  pre¬ 
dictions  for  SOD  and  sediment  release  of  NH4,  NO3,  and  DIP  are  shown 
in  Figure  15.  The  results  show  that  increased  deep  burial  velocities  pro¬ 
duced  lower  SOD  and  sediment-water  fluxes  of  NH4,  NO3,  and  DIP.  The 
sediment  release  of  nitrogen  and  phosphorous  was  reduced  at  the  sites 
where  the  deep  burial  rate  was  increased.  However,  the  magnitude  of  the 
burial  velocity  influenced  by  the  depositional  flux  of  solids  from  the  water 
column  and  the  compaction  are  not  included  in  the  current  model. 


Figure  15.  Sensitivities  of  deep  burial  velocities  on  the  sediment-water  fluxes:  (a) 
SOD,  (b)  NH4,  (c)  N03,  and  (d)  DIP. 
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4  Conclusions 

A  benthic  sediment  diagenesis  module  has  been  developed  and  included  in 
the  NSM  modules  as  an  option.  This  new  feature  required  additional  mod¬ 
el  inputs.  When  this  feature  is  activated,  SOD  and  water-sediment  fluxes 
are  internally  computed  by  the  sediment  diagenesis  module  for  coupling 
with  water  column  kinetics.  Three  major  processes  included  in  the  sedi¬ 
ment  diagenesis  module  are:  1)  deposition  of  particulate  organic  matter 
from  the  water  column  to  the  sediments,  2)  mineralization  (or  diagenesis) 
of  the  particulate  organic  matter,  and  3)  reactions  and  transfers  of  the  re¬ 
action  products. 

Testing  and  validation  studies  of  the  NSMII-benthic  sediment  diagenesis 
module  were  performed  by  running  the  same  site  conditions  with  the  Sed- 
Flux  and  ICM  for  the  Chesapeake  Bay  dataset.  In  general,  the  model  verifi¬ 
cations  with  respect  to  temporal  variations  of  sediment  organic  and 
inorganic  carbon,  nitrogen,  and  phosphorus  concentrations  are  quite  satis¬ 
factory  compared  with  results  produced  from  the  SedFlux  and  ICM.  The 
NSMII  was  able  to  reproduce  SOD  and  sediment-water  nutrient  fluxes. 
Moreover,  the  sediment  diagenesis  module  was  applicable  to  a  wide  range 
of  sedimentation  conditions  and  sediment  organic  contents.  There  was 
good  agreement  between  the  NSMII  predictions  and  results  produced 
from  the  SedFlux  and  ICM  for  SOD,  sediment-water  fluxes  of  NH4,  NO3, 
DIP,  H2S  and  CH4.  Special  test  cases  and  sensitivity  analysis  of  the  NSMII 
to  the  sediment  initial  conditions,  organic  matter  depositional  rates,  deep 
burial  velocities,  and  particle  mixing  coefficient  highlighted  process  inter¬ 
actions  and  synergistic  effects.  It  is  concluded  that  the  NSMII-benthic  sed¬ 
iment  diagenesis  module  will  be  most  useful  for  applications  involving 
sediment- water  interactions  to  evaluate  the  effects  of  benthic  sediments 
on  water  quality  concentrations  in  aquatic  environments. 

There  are  some  limitations  of  the  NSMII-benthic  sediment  diagenesis 
module  that  could  be  improved  in  the  future.  The  benthic  sediment  layer  is 
assumed  to  have  constant  properties,  including  the  thickness,  volume,  po¬ 
rosity,  bulk  density,  and  solids  mass.  The  NSMII  does  not  simulate  or  cou¬ 
ple  with  sediment  erosion  and  resuspension  processes.  Particle  deposition 
and  erosion  and  the  resulting  effects  on  bed  thickness  are  not  taken  into 
account.  Coupling  the  NSM  module  with  a  dynamic  benthic  sediment 
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model  within  the  hydrologic  and  hydraulic  models  could  overcome  some 
of  the  weaknesses. 
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